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Sudden transition from naked atom decay to dressed atom decay 
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The studies on quantum open system play key roles not only in fundamental problems in quantum 
mechanics but also in quantum computing and information processes. Here we propose a scheme 
to use a one dimensional coupling cavity array (CCA) as an artificial electromagnetic environment 
of a two-level atom. For a finite length of CCA, we find that after a turning time the population of 
excited state deviates suddenly from the exponential decay. We show that physically this phenomena 
corresponds to a transition from a naked atom decay to a dressed state decay. We hope that our 
finding will promote the studies on quantum system with a finite size environment. 

PACS numbers: 42.50.Pq, 32.80.Qk, 03.65.Yz 


I. INTRODUCTION 

In nature there does not exist a completely closed 
quantum system. Any quantum system must interact 
with environment, and the environment in most cases 
is very complex and often not precisely specified, which 
makes it extremely difficult to study such an open sys¬ 
tem. One obvious reason for us to study an open system 
is that the knowledge on the system’s state can not be 
obtained without the interaction between the measure¬ 
ment apparatus and the system. In addition, to control 
a quantum system to implement some quantum comput¬ 
ing and information processes, we have to fight against 
the decoherence and dissipation induced by the environ¬ 
ment L lj. 

Our aim is to study the effect of an environment 
with finite size. We model an environment with a one¬ 
dimensional coupling cavity array, which may be real¬ 
ized in experiments by coupled superconducting trans¬ 
mission line resonators or defect resonators in photonic 
crystals M- In literature, a one-dimensional coupling 
cavity array with a two-level atom in is first introduced to 
study single photon scattering problem j6|. Then many 
variants of this model are studied in recent years, for ex¬ 
ample, a super cavity can be formed by embedding two 
atoms to control the tunneling or directly decreasing two 
tunneling strengths [J [§] • In addition, when putting an 
excited atom in the middle cavity of the CCA, the atom 
finally will evolve into two bound states and illustrate the 
oscillation behavior 

Here we study the spontaneous decay of a two-level 
atom located in one node of a one dimensional (ID) cou¬ 
pling cavity array (CCA). For the atom, the ID CCA can 
be regarded as a ID electromagnetic environment. Be¬ 
cause the length of the cavity array can be controlled in 
experiments, we have the opportunity to study the effect 
of an environment with finite size. For an environment 
with finite size, we may expect the spontaneous decay of 


the atom will not be exponential one EMI- We find 
that there is a sudden deviation of the exponential decay 
at some turning time. We theoretically explain this phe¬ 
nomena as the transition from a naked atom decay to a 
dressed atom decay, and numerical results confirm this 
underlying mechanism. 

The remainder of the paper is organized as follows. In 
Sec. 2, the theoretical model of a two-level atom in a ID 
CCA is introduced as a microscopic model of an open 
quantum system. In Sec. 3, we analyze the condition 
for the exponential decay in our system, and gives the 
analytical result on the decay rate. In Sec. 4, we give 
the numerical results on the atom decay for the CCA 
with finite length, and find that the deviation of the ex¬ 
ponential decay is a sudden transition at some turning 
time. We give a physical explanation why this happens 
and predict the turning time. In Sec. 5, after the turning 
time, we gives the results based on the dressed atom de¬ 
cay, which is confirmed by the numerical results. Finally 
we will give some discussions and a brief summary. 

II. MODEL: A TWO-LEVEL ATOM IN ID CCA 

We consider a two-level atom located in the n-th cavity 
of a ID CCA with length N as shown in Fig. []] where 
the atom is our quantum system and the CCA is its ID 
electromagnetic environment. 
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FIG. 1. (Color online). Schematic configuration of our sys¬ 
tem. A two-level atom (filled blue circle) is in the n-th cavity 
of the CCA. 

Under the rotating wave approximation the Hamilto¬ 
nian of the system becomes 
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U = Ho + Th + Ha 


(1) 
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where the Hamiltonian of the CCA Hq is 
N N 

U 0 =UcY, a |fflj - y (gtfflj-i + h.c.j (2) 

3=1 3=1 

with uj c the frequency of each single-mode cavity, 77 the 
hopping strength, N the number of cavities in the CCA, 
and aj (rij ) the photon creation (annihilation) operator 
for the j-th single-mode cavity. Hereafter, we set 77 = 1 
and w c = 0; The Hamiltonian for the free two-level atom 
7i a is 

7-La = w 0 |e)(e|, (3) 

where w Q is the energy of the excited state, \e) (|<?)) is 
the excited (ground) state of the atom. The interaction 
between the atom the cavity 7~Li is 

7ii =9 (ala- + a n cr+) , (4) 

where g is the coupling the atom and the cavity, and 
u _ (cr + ) is the atomic lowering (raising) operator. 

Notice that the total photon number n p = (CjLi a ] a j 
in H 0 is conserved for the CCA. In the subspace of n p = 
1, the eigenstate and eigenvalue of the CCA is 


Due to the excitation number is conserved in our sys¬ 
tem, the state at any time t can be expanded as 

N 

l'f'(t)) = [a{t)a + + ^2p k (t)b]]\vac)\g), (8) 

fc =1 

where a{t) and /3k(t) represent the excitation probability 
amplitude for the atom state and the k- th mode at time 
t respectively. 

Through the Schrodinger equation 

iA|*( t ))= W |tt(t)>, (9) 

we get the following set of equations about the coefficients 

N 

ia(t) = - iui a a(t ) - i^gkPkjt), ( 10 ) 

fc =1 

iPk(t) = -iwfc/3fc(<) - igkOt{t). ( 11 ) 

Making the transformation of parameters 
5(f) = a{t)e iuiat 

J k {t) = pkity^ ( 12 ) 

A k =U> k - IjJa 


/ 2 \ 

Wk) = J sm{j9 k )a]\vac) = b\\vac), (5) 

3=1 

uj k = -2cos0 fc , (6) 

where k is integer with 1 < k < N, 9 k = Jfbi-i b k ( 6 ^) the 
the photon creation (annihilation) operator for the fc-th 
mode and |uac) is the vacuum state of the cavity array. 

Using Eqs. ([5]) and d 6 j) , we can rewrite the Hamiltonian 
as 

N N 

Ti = ^cjkbjbk + ^2g k (b{cr- + h.c.j +w a |e)(e| (7) 
k =1 fc =1 


with g k = s ia( n ®k)9 being the coupling between 

the atom and the fc-th mode and. Notice that this is a 
typical spin-boson model to study quantum open system, 
and here the energy spectrum oj k of the bosonic environ¬ 
ment and the coupling strength g k are specified. 


III. EXPONENTIAL DECAY OF THE NAKED 
ATOM STATE 


Now we study the spontaneous decay of the two-level 
atom in our system. In other words, when the atom is 
initially prepared in the excited state and the CCA stays 
in its vacuum state, what is dynamics of the excitation 
probability of the atom? In this section, we will focus on 
studying the conditions for the exponential decay of the 
excited atom, and theoretically predict the decay rate. 


we simplify the equations of parameters as 

N 

a(t) =g k p k (t)e~ zAkt (13) 

k= 1 

/3fc(f) = -ig k a{t)e lAkt (14) 

Integrate Eq. (lldl) then substitute it into Eq. m , we 
obtain an equation only about a(t) 


N 


= - Z) 


9k 


dt'a(t')e 




(15) 


fc =1 


The exponential decay can be achieved when the near 
resonant modes are closely spaced in frequency and share 
almost the same coupling strength with the atom, mean¬ 
while the other modes far off resonance are of little im¬ 
portance in the decay. To meet up with the requirements, 
we can use a cavity array with length long enough and 
put the atom in a selected location. Then the decay rate 
can be analytically derived. In our model, it is 


r = 


4 g 2 


(16) 


with ujo = 2 cos k 0 = u} a the energy of the resonant mode 
(see Appendix for the detailed calculation). 

We now give some numerical results to support the the¬ 
oretical analysis above. Figure[5])a) shows the modes con¬ 
dition around the resonant frequency. With the atom set 
in the 1984-th cavity (the antinode of the resonant mode), 
the above mentioned requirements are satisfied. Further¬ 
more, in a short time the behavior of the atom state is a 
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FIG. 2. (Color online), (a) The 50-th to 60-th modes of CCA 
with N = 2001. The yellow solid line stands for the mode (55- 
th) resonant with the atom while the other blue solid lines 
are for the modes nearby, (b) Time dependent atom state 
excitation for various lengths of CCA. The solid lines are the 
exact numerical results as for red N = 1001, n = 992, black 
N = 1501, n = 1488 and green N = 2001, n = 1984. The 
corresponding dash lines are exponential decay lines with the 
decay rate determined by Eq. (11611 . Here, g = 0.0015 and the 
atom is resonant with the 55-th mode in all situations. 



FIG. 3. (Color online), (a) Black dash line is the atom exci¬ 
tation and solid lines stand for various modes, (b) Only the 
near-resonance modes, for the solid lines each color represents 
the same mode as in (a), the yellow dashed line stands for the 
results calculate from Eq. m- To simplify the figure, we 
only show some blue tuned near-resonance modes while not 
including the red tuned ones, their behavior actually are very 
much alike. 


truly well-defined exponential decay and with its decay 
rate well predicated by Eq. (fl6l) (see Figure [U(b)). 

But when time passes a specific point, the evolution 
dramatically deviates from the exponential line, and the 
change can be postponed by lengthening the cavity array. 
Next, we will explain the appearance of the turning point, 
and show that after the point a new physical process 
starts. 


IV. DEVIATING FROM THE EXPONENTIAL 
DECAY 


with A k = t2/ 9 } +uj , . 

Eq. (fl9l) can perfectly describe the evolution of the 
modes during the time when the atomic state obeys the 
exponential decay as shown in Fig[3](b). For A*,, it sat¬ 
isfies the following equation A k « (k — fco)Ai (see Ap¬ 
pendix for proof), here we define Ai = oj ko +i — uj ko ■ To 
determine the minimum of M kl we calculate its derivative 

M k (t) = A k e- Tt / 2 {-Ye- Tt ' 2 + r cos A k t + 2 AfcsinAfc). 

( 20 ) 

For A fc T, only the third term in the right-hand side 
of Eq. (12011 is predominant. Then minimums appear near 


The key point for understanding the deviation is to 
clarify the behavior of the modes around the resonant 
frequency. We pick out the N = 2001 condition for exam¬ 
ple, Fig[3](a) reveals the evolution of the modes alongside 
with the atom state. Around the turning point, except 
the resonant mode other near-resonance modes almost 
all together oscillate to near zero. Fig [51(b) gives a more 
clear view about this pattern. This is the main reason to 
cause the deviation since the resonant mode becomes pre¬ 
dominate over other near resonant modes at the turning 
point, which violates the requirements for the exponen¬ 
tial decay of a naked atom. 

Now we turn to give analytical analysis about the ap¬ 
pearance of this pattern. As the atom state exponentially 
decay with time, its excitation probability amplitude is 

a(t) = e _rt / 2 , (17) 


Then the k -th mode’s excitation probability amplitude 
can be determined through Eq. m 


hit) = - ig k f dt'e^- T / 2 V 
Jo 

(iA fc -r/2 )t _ -q 




(18) 


iu k - T/2 

while its probability is 


M k (t) = \/3 k (t)\ 2 = A k {e rt — 2e rt/2 cos A k t + 1), (19) 


27 tI 27r/ 

A k \k-k 0 \A 1 


( 21 ) 


as l is positive integer. 

With Eq. (1211) we can see when the excitation of the 
nearest mode (mode ko ± 1 ) reach the first minimum, 
the other near-resonance modes also reach its minimum 
(though they have already oscillated several times). Then 
the time of the turning point t c is totally dependent on 
the energy difference between the resonant and nearest 
near-resonance mode 


t c = 


2tt 

aT 


( 22 ) 


And this is well consistent with the numerical results. 
When lengthening the cavity array, Ai decreases and 
leads to the postpone the turning point as we mentioned 
before. For the extreme case, a ID CCA with infinite 
length Ai —0 leads to t c —> oo. So this deviation phe¬ 
nomenon is actually due to the finite size of the environ¬ 
ment. 


V. DECAY OF THE DRESSED ATOM 

Around the turning point, the effect of the reso¬ 
nant mode will be predominant since the other near¬ 
resonance modes are almost empty without any photon. 
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We conclude that the atom and the resonant mode will 
be str ong ly coupled to form a subsystem, the dressed 
atom So for the next stage of time, it is the decay 
of the dressed atom, a physical process totally different 
with before. 

We can introduce a general master equation to describe 
the evolution of the dressed atom 0 

d r 

P = -i[H d ,p\- ^ (<J + cr~p d + pcr + a~)+Ta~pa + (23) 
with 

Hd = cr z + iOob\b r + g r {b\a~ + h.c.) (24) 


the Hamiltonian of the dressed atom, g r {= y n+\ 9) 
the coupling strength between the atom and the resonant 
mode, p the density matrix of the dressed atom and T the 
decay rate of the subsystem (same as the the spontaneous 
emission rate of the atom). 

Under the secular approximation, we can project the 
master equation over the dressed states 

d r 

jjPii = ^ Pa (25) 

d r 

~^P12 = (-2 ig r - -)pi2 (26) 

with i = 1,2, Pij = (i\p\j), |1) = -±=(\vac)\e) + |lfc 0 >|s» 
and |2) = ^(|uac)|e) — |lfc 0 )|<7)) the dressed states. 
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FIG. 4. (Color online). Atom excitation as a function of 
time. Solid lines represent exact numerical results with differ¬ 
ent color stands for different cavity array length. The dotted 
lines are obtained through Eq. (1281 . All the parameters are 
the same as before. 

For the atom excitation 

Pee(t) = 2 [PllW + P22(t) + Pl2 (t) + P2l(t)\ (27) 


Substitute the results obtained through Eq. (l25l) and 
Eq. (l26l) into Eq. (l27l) . the evolution of the atom state is 
determined as 

Pee{t) = 2 e C \pil{ic) + P 22 (t c ) _ 

+ Pi2(t c )e ~ 219 '"<*-*4 + P 2 i(t c )e 2l9 ^ t -^] 

with p ee (t) = {e\{vac\p(t)\vac)\e). 

Fig 0] shows that after turning point, the dressed atom 
decay theory can describe the system evolution quite ac¬ 
curately. But after the time point when the atom state 
is empty, the dressed atom system no longer exists, the 
evolution becomes quite chaotic and unpredictable. 

In conclusion, we study the atom decay within the ID 
CCA frame. By selecting appropriate parameters, we 
achieve the well-known exponential decay of the atom 
state, and theoretically give the decay rate. Meanwhile 
the exponential decay only last a limited time in this 
model, the evolution deviates after a time point we name 
t c . This is due to the near-resonance modes almost be¬ 
come empty together around t c , and we theoretically 
prove the existence of t c in our model. After t c , the 
resonant mode is strongly coupled with the atom to form 
the dressed atom subsystem, and the dressed atom de¬ 
cay theory can describe the evolution quite accurately. 
Finally when the atom evolves into its ground state, the 
dressed atom subsystem no longer exists, and the evolu¬ 
tion after is quite chaotic and beyond our understanding 
now. We hope our work can enlighten the study of the 
transition between the atom state exponential decay and 
dressed atom decay, and help to have a better under¬ 
standing about the dynamics of quantum open system 
with a finite size environment. 
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APPENDIX: HOW MANY MODES ARE 
ENOUGH? 

Throughout the paper, we emphases that only modes 
around the resonant frequency are important in the sys¬ 
tem evolution while the far off-resonance modes can be ig¬ 
nored (this is also the key point in the Weisskopf-Wigner 
approximation). Then the question appears immediately, 
how to define near-resonance and far off-resonance or in 
another way say how many modes is actually enough? 
We give some numerical results to show intuitively that 
the number is quite small. 

For the cavity array with N = 2001, Fig [5] reveals that 
the system can be described quite accurately using only 
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5 modes and perfectly with 11 modes. Compared with 
the total 2001 modes in the system, this is a very small 
number. 



Adopting the Weisskopf-Wigner approximation mm, 
we get 


pOO p 2 r\ 2 

5(t) = —a(t) / dr / dw— 7 =L=e < ("»-“) T (32) 
J o J-2 7 TV 4 — 

with r = t — t' . Since 


/•°° i 

/ dre- i( “-“ o)T = 7r6(td - u 0 ) - »P(-) (33) 

Jo W — Wo 

where P represents the Cauchy principal part which leads 
to the well-known lamb shift. Since the lamb shift has 
nothing to do with decay, we neglect it here and obtain 

a(t) = e ~** (34) 

with r the atom spontaneous emission rate 


FIG. 5. (Color online). Time dependent atom excitation. The 
green solid line is the exact results including all the modes, 
the dotted blue line including the 50-th to 60-th modes and 
the red dashed line including the 53-th to 57-th modes. Here 
the other parameters are the same as before. 


r = 


4 g 2 3 

\/ 4 ^ w o 


(35) 


APPENDIX: CALCULATION OF THE DECAY 
RATE 


APPENDIX: PROOF OF A k « {k - fc 0 )Ai 


Prom Eq. (ED and Eq. ©, we have 


As mentioned above, for the modes around the reso¬ 
nant frequency the coupling strength g k « g r = \J n+i 9- 
then we can rewrite Eq. m as 

o n 2 N pt 

«(*) = E / dt'5(t')e- lAfe(t - 4} . (29) 

+ 1 k = i Jo 

After replacing the summation over k with an integration 
as d k = we get 

S(t) = -^ [" dk f dt'ait’y- 1 ^-^. (30) 

7T Jo Jo 


. , kir k 0 TT 

Afc = 2 (cos ——7 — cos ■ 


N + 1 


N + l' 


(36) 


Using the sum and difference formulas of the trigonomet¬ 
ric functions, we rewrite it as 


Afc 


—4 sin 


rmr 

2(N + 1) 


sin 


2ko + m 
2(N + 1) ’ 


(37) 


where m = k — ko- Since we only consider the modes 
around the resonant frequency, m <C N and m <C k$. 
Then Eq. (1371) can be simplified as 


By replacing the integral variable k with frequency to 

since u = 2 cos k, we have Afc ~ ~ 4 2 (a + 1 ) Sin 2{N + +1) ” ^ ~ fc ° )Al ' (38) 

5(t) = -H£! [* dt' [ du , 1 

7T Jo J —2 v4-w 2 
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